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Phase diagram of strongly interacting matter is discussed within the exactly solvable statistical 
model of the quark-gluon bags. The model predicts two phases of matter: the hadron gas at a low 
temperature T and baryonic chemical potential /is, and the quark-gluon gas at a high T and/or 
fiB ■ The nature of the phase transition depends on a form of the bag mass- volume spectrum (its 
pre-exponential factor), which is expected to change with the Hb/T ratio. It is therefore likely that 
the line of the I s * order transition at a high (j,b/T ratio is followed by the line of the 2 nd order 
phase transition at an intermediate hb/T, and then by the lines of "higher order transitions" at a 
low [LB IT. 
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Investigation of the properties of strongly interacting matter at a high energy density is one of the most important 
. subjects of the contemporary physics. In particular, the hypothesis that at high energy densities the matter is in 
' the form of quark-gluon plasma (QGP) 1] rather than a gas of hadrons (HG) motivated a first stage of the broad 
I/"") , experimental program of study of ultra-relativistic nucleus-nucleus collisions Q ■ Over the last 20 years rich data were 

■ collected by experiments located at Alternating Gradient Synchrotron (AGS) and Relativistic Heavy Ion Collider 
\ (RHIC) in Brookhaven National Laboratory, USA and at Super Proton Synchrotron (SPS) in CERN, Switzerland. 

■ The results indicate that the properties of the created matter change rapidly in the region of the low SPS energies Q 
(V^nn ~ 6-12 GeV). The observed effects confirm predictions for the transition from hadron gas to QGP Q and thus 
indicate that in fact a new form of strongly interacting matter exists in nature at sufficiently high energy densities. 

What are the properties of the transition between the two phases of strongly interacting matter? This question 
motivates the second stage of the investigation of nucleus-nucleus collisions. A new experimental heavy ion program 
O ' (V^nn ~ 5-20 GeV) at the CERN SPS is in preparation ||| and a possible study at the BNL RHIC is under discussion. 

Based on the numerous examples of the well-known substances a conjecture was formulated (see review paper 
Ph ■ @ and references there in) that the transition from HG to QGP is a I s * order phase transition at low values of a 
temperature T and a high baryo-chemical potential \xq and it is a rapid cross-over at a high T and a low hb- The 
end point of the I s * order phase transition line is expected to be the 2 nd order critical point. This hypothesis seems 
/\ ' to be supported by QCD-based qualitative considerations [|| and first semi-quantitative lattice QCD calculations 0. 
^ . The properties of the deconfinement phase transition are, however, far from being well established. This stimulates 
our study of the transition domain within the statistical model of quark-gluon bags. 

The bag model [8( was invented in order to describe the hadron spectrum, the hadron masses and their proper 
volumes. This model is also successfully used for a description of the deconfined phase (see e.g. Ref. |!j). Thus the 
model suggests a possibility for a unified treatment of both phases. Any configuration of the system and, therefore, 
each term in the system partition function, can be regarded as a many-bag state both below and above the transition 
domain. The statistical model of quark-gluon bags discussed in this paper contains in itself two well-known models of 
deconfined and confined phases: the bag model of the QGP [fj and the hadron gas models [T^j. They are surprisingly 
successful in describing the bulk properties of hadron production in high energy nuclear collisions and thus one may 
hope that the model presented here may reflect the basic features of nature in the transition domain. 

The paper is organized as follows. The important properties of the statistical bootstrap model with the van der 
Waals repulsion are summarized in Sect. II. In Sect. Ill the statistical model of the quark-gluon bags is presented. 
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Within this model properties of the transition region for /j,b = are studied in Sect. IV and the analysis is extended 
to the complete T — ps plane in Sect. V. The paper is closed by the summary and outlook given in Sect. VI. 



II. STATISTICAL BOOTSTRAP MODEL AND VAN DER WAALS REPULSION 



The grand canonical partition function for an ideal Boltzmann gas of particles of mass m and a number of internal 
degrees of freedom (a degeneracy factor) <?, in a volume V and at a temperature T is given by: 



N=0 j=l y ' 
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where 



(f>(T, m) 



1 

2^2 



k 2 dk exp 



(k 2 



? 2U/2 



T 



m 2 T 



and Ki is the modified Besscl function. The function 4>(T, to) is equal to the particle number density: 

_ N(y,T) 

H T ) = — y — = 9<p{T,m) . 

The ideal gas pressure and energy density can be derived from Eq. (yi as: 
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One can easily generalize the ideal gas formulation £[J to the mixture of particles with masses mi,...,m„ and 
degeneracy factors g% , . .. , g n : 



Z(V,T) 



E ••• E 

AT 1= N n =0 



N-l\ 



A'i 



[Vg n <p(T, m n )} 
NJ. 



N n 



exp 



V Es^^V' 

3 = 1 



(5) 



The sum over different particle species j can be extended to infinity. It is convenient to introduce the mass spectrum 
density p(m), so that p(m)dm gives the number of different particle mass states in the interval [m,m + dm], i.e. 
YlJLi 9j--- = Jo°° dm...p(m). In this case the pressure and the energy density are given by: 



P(T) 
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Jo 

Eqs. P were introduced by Hagedorn |Tlj for the mass spectrum increasing exponentially for m — > oo: 

p(m) ~ C m" a exp(6m) , (7) 

where a, b and C are model parameters. This form of the spectrum Q was further derived from the statistical boot- 
strap model It can be shown, that within this model the temperature Th = 1/6 (the "Hagedorn temperature") 
is the maximum temperature of the matter. The behavior of thermodynamical functions 10 with the mass spectrum 
(|7|) depends crucially on the parameter a. In particular, in the limit T — > Th the pressure and the energy density 
approach: 



p , e -> oo , 
p — * const , e 
p , e — > const 



for a < - ; 
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(8) 
(9) 
(10) 



Up to here all particles including those with to — > oo were treated as point-like objects. Clearly this is an unrealistic 
feature of the statistical bootstrap model. It can be overcome by introduction of hadron proper volumes which 
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simultaneously mimic the repulsive interactions between hadrons. The van der Waals excluded volume procedure 0] 
can be applied for this purpose 0] (other approaches were discussed in Refs. 0,0). The volume V of the ideal 
gas is substituted by the "available volume" V — v a N, where v Q is a parameter which corresponds to a particle 
proper volume. The partition function then reads: 

*wt)- f » f -y"i% (v-a). an 

N=0 

The pressure of the van der Waals gas can be calculated from the partition function l|llfl by use of its Laplace 
transform. This procedure is necessa ry b ecause the "available volume" , V — v a N, depends on the particle number N. 
The Laplace transform of JTTJ is fjUllfij: 

Z{s,T) ee [°° dVexp(-sV) Z(V,T) = V f° dV (y _ VoN f exp( - sV ) 



1 oo 

jv=o 
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exp(-« oS ) g<f>(T,m)} 1 . (12) 



In the thermodynamic limit, V — » oo, the partition function behaves as Z(V,T) ~ exp[pV/T]. An exponentially 
increasing Z(V,T) generates the farthest-right singularity s* = p/T of the function Z(s,T) in variable s. This is 
because the integral over V in Eq. I|12|) diverges at its upper limit for s < p/T. Consequently, the pressure can be 
expressed as 

P (T) = T lira h^lll = T s*(T) , (13) 

and the farthest-right singularity s* of Z(s,T) (|12|l can be calculated from the transcendental equation 0,0]: 

s*(T) = exp[-t; s*(T)] g4>(T,m) . (14) 

Note that the singularity s* is not related to phase transitions in the system. Such a singularity exists for any 
statistical system. For example, for the ideal gas (v — in Eq. I|14[l) s* = g<p(T, m) and thus from Eq. (|13fl one gets 
p = Tg<f>(T, m) which corresponds to the ideal gas equation of state (@J. 



III. GAS OF QUARK-GLUON BAGS 



The van der Waals gas consisting of n hadronic species, which are called bags in what follows, is considered in this 
section. Its partition function reads: 
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(15) 



where (mi, «i), (m„, v n ) are the masses and volumes of the bags. The Laplace transformation of Eq. Ijl5(l gives 



Z(s,T) = 



- ^2exp(-vjs) g (j){T,mj) 

3 = 1 



(16) 



As long as the number of bags, n, is finite, the only possible singularity of Z(s, T) (jlfijl is its pole. However, in the case 
of an infinite number of bags the second singularity of Z(s, T) may appear. This case is discussed in what follows. 

Introducing the bag mass- volume spectrum, p(m,v), so that p(m,v)dmdv gives the number of bag states in the 
mass- volume interval [m, v; m + dm, v + dv], the sum over different bag states in definition of Z(V, T) can be replaced 
by the integral, J^JLi 9j--- = Jq° dmdv...p(m,v). Then, the Laplace transform of Z(V,T) reads [l3j |: 



Z(a,T) 



dVexp(-sV) Z(V,T) = [ 
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where 

f(T,s) = I dmdv p(m,v) exp(— vs) 4>(T,m) . (18) 



The pressure is again given by the farthest-right singularity: p{T) = Ts*(T). One evident singular point of Z(s,T) 
(TPTIl is the pole singularity, sh{T): 

s H {T) = f(T,s H (T)) . (19) 

As mentioned above this is the only singularity of Z(s,T) if one restricts the mass-volume bag spectrum to a finite 
number of states. For an infinite number of mass- volume states the second singular point of Z(s,T) l|17|). sq(T), can 
emerge, which is due to a possible singularity of the function f(T, s) <|18|) itself. The system pressure takes then the 
form: 

P (T) = Ts*(T) = T-max{s H (T),s Q (T)} , (20) 

and thus the farthest-right singularity s*(T) of Z(s, T) 11170 can be either the pole singularity sh(T) (11911 or the sq(T) 
singularity of the function f(T, s) (|18|l itself. The mathematical mechanism for possible phase transition (PT) in our 
model is the "collision" of the two singularities, i.e. sh{T) = sq(T) at PT temperature T = Tc (see Fig. 1 . In 
physical terms this can be interpreted as the existence of two phases of matter, namely, the hadron gas with the 
pressure, pn — Tsn(T), and the quark gluon plasma with the pressure pq = Tsq(T). At a given temperature T 
the system prefers to stay in a phase with the higher pressure. The pressures of both phases are equal at the PT 
temperature Tc- 

An important feature of this modeling of the phase transition should be stressed here. The transition, and thus the 
occurrence of the two phases of matter, appears as a direct consequence of the postulated general partition function 
(a single equation of state). Further on, the properties of the transition, e.g. its location and order, follow from the 
partition function and are not assumed. This can be confronted with the well-known phenomenological construction 
of the phase transition, in which the existence of the two different phases of matter and the nature of the transition 
between them are postulated. 

The crucial ingredient of the model presented here which defines the presence, location and the order of the PT is 
the form of the mass- volume spectrum of bags p{m, v). In the region where both m and v are large it can be described 
within the bag model [8j . In the simplest case of a bag filled with the non- interacting massless quarks and gluons one 
finds IH : 



p(m,v) ~ C v 1 {m — Bv) s exp 



4 



r y* ,,V4 
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(21) 
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where C, 7, 5 and B (the so-called bag constant, B w 400 MeV/fm 3 Q) are the model parameters and 

■m-s (* + !**) =fiK 2 ' 8 + r 2 ' 2 - 3 - 3 ) =55? < 22 > 

is the Stefan-Boltzmann constant counting gluons (spin, color) and (anti-)quarks (spin, color and u, d, s-flavor) 
degrees of freedom. This is the asymptotic expression assumed to be valid for a sufficiently large volume and mass of 
a bag, v > Vq and m > Bv + Mo- The validity limits can be estimated to be Vo « 1 fm 3 and Mo « 2 GeV [20|. The 
mass- volume spectrum function: 

n 

p H (m,v) = gj 5(m - rrij) S(v - Vj) (23) 

should be added to p(m, v) in order to reproduce the known low-lying hadron states located at v < Vo and m < 
BVo + Mq- The mass spectra of the resonances are described by the Breit-Wigner functions. Consequently, a general 
form of f(T, s) (JTHJ reads: 



n r DO r OC 

f(T,s) = / g (T, s) + f Q (T, s) = gj exp(-Vjs) cj>(T, mj) + dv dm 

j =1 Jv o J M„+Bv 



p(m,v) exjp(—sv)(j>(T,m) , (24) 



intermediate energy nucleus-nucleus collisions 
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where p(m, v) is given by Eq. J5TJ. 

The behavior of Jq(T,s) is discussed in the following. The integral over mass in Eq. (|24|) can be calcu- 
lated by the steepest decent estimate. Using the asymptotic expansion of the /^-function one finds 4>(T, m) ~ 
(mT/27i") 3//2 exp(— m/T) for m^$> T. A factor exponential in to in the last term of Eq. i|24[l is given by: 



exp 



exp[[/(m)] 
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Presenting [/(to) as 
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one finds 



(—) 

\dm 2 I 



SQ{T) 



(to - to ) 2 = vsq(T) 



(to — Too) 2 
8<jq vT 5 



B 
T 



f Q (T,s) ~ u(T) 
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(25) 
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Vi, 



where u(T) = C%- 1 a S Q +1/2 T 4+4S (<tqT 4 + B) 3/2 . The function f Q (T,s) has the singular point s = s Q because 
for s < sq the integral over dv diverges at its upper limit. 

The first term of Eq. 124fl . fjj, represents the contribution of a finite number of low- lying hadron states. This 
function has no s-singularities at any temperature T. The integration over the region m, v — * oo generates the 
singularity sq (T) (|28|l of the fg (T, s) function (|29|l . As follows from the discussion below this singularity should be 
associated with the QGP phase. 

By construction the function f(T, s) (|18|l is a positive one, so that sjy(T) l|19|) is also positive. On the other hand, 
it can be seen from Eq. (|2"S|) that at a low T sq(T) < 0. Therefore, s H > sq at small T, and according to Eq. (|2"Ut it 
follows: 



/'I T ) : T,„ . :(/) 2 ^ 



(30) 



The system of the quark-gluon bags is in a hadron phase. 

If two singularities "collide", T = Tc and sh{Tc) = sq(Tc), the singularity sq(T) can become the farthest-right 
singularity at T > Tc- In Fig. ^ the dependence of the function f(T,s) on s and its singularities are sketched for 
T\ < T2 = Tc < T3. The thermodynamical functions defined by the singularity sq(T) are: 



P(T) 



= T sn = — T 4 
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ds 



e(T) = T'=g- = <r Q T 
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(31) 



and thus they describe the QGP phase. 

The existence and the order of the phase transition depend on the values of the parameters of the model. From 
Eq. (|2~9"|l follows that a PT exists, i.e. s* — sq at high T, provided 7 + <5 < — 3 (otherwise f(T, sq) = 00, and s# > sq 
for all T, for illustration see Fig.|2 left). 

In addition it is required that f(T, sq) < sq(T) at T — > 00, (otherwise s* = sh > sq for all T, for illustration see 
Fig. [21 right). For 7 + 6 < -3 one finds f Q {T, sq) cx T 10+4S at T -► 00. On the other hand, s Q {T) cx T 3 at a high T 
and therefore 8 < —7/4. Consequently the general conditions for the existence of any phase transition in the model 
are: 



7 < 



S < - 



7 



(32) 



The parameters 7 and S can be calculated within the model provided constraints defining allowed quark-gluon bags are 
given |l9l |20| . Assuming massless non-interacting quarks and gluons with the condition of the total three-momentum 
vector equal to zero one gets fl9| : 



7 = 



A 7 

6 = 4' 



(33) 



G 




FIG. 1: The dependence of f(T,s) on s for three different temperatures: T\ < T2 = Tc < T3 (solid lines). The pole 
singularity sh and the singularity sq are denoted by circles and crosses, respectively. A PT corresponds to the "collision" of 
two singularities sh = sq at the temperature Tc- 




FIG. 2: The dependence of f(T,s) on s for two different temperatures, Ti < T2 (solid lines). The figures demonstrate the 
absence of the PT in the model if f(T, sq) — 00 (left), or if f(T, sq) > sq at all T (right). 



and consequently there is no phase transition in the model. For a graphical illustration of this case see Fig- El left. An 
additional selection of only colorless bags yields |20|: 

9 19 

7 = " J , 5 = ~ T ' (34) 

so that the conditions ll'>2[| are fulfilled and the model leads to the phase transition, see Fig. 1 for illustration. Based 
on these two simple examples one concludes that the parameters 7 and d depend significantly on the set of selected 
constraints. These constraints in statistical mechanics are usually ignored as the related to them pre-exponential 
factors do not contribute to the equation of state in the thermodynamic limit. Because of this possible sets of physical 
constraints are not well defined and their consequences were not studied. Therefore in what follows 7 and i5 are 
considered as free model parameters. 
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IV. 1 st , 2 nd AND HIGHER ORDER PHASE TRANSITIONS 

In this section the order of the PT in the system of quark-gluon bags is discussed. 
The 1 st order PT takes place at T = T c provided: 



s Q {T c ) = s H (T c ) , 



> 



T=T C 



(ds H \ 
\dT ) T= 



(35) 



T=T C 



Thus the energy density (e = T 2 ds/dT) has discontinuity (latent heat) at the 1 st order PT. Its dependence on T is 
shown in Fig. [31 left. In calculating dsu/dT it is important to note that the function sh{T) (|19fl is only defined for 
T<T C , i.e. for s„(T) > s Q (T). 






FIG. 3: The dependence of the energy density e on temperature for the 1 st (left) and 2 nd (right) order PT. The energy density 
e(T) has a discontinuity at T = T C for the 1 st order PT whereas for the 2 nd order PT e(T) is a continuous function of T, but 
de/dT has a discontinuity at T = Tc- 



The 2 nd order PT takes place at T = T c provided 
s q( t c) = s H (T c ) , 



ds Q 
dT 



T=T C 



ds H 
dT 



T=T C 



d 2 s H 
d 2 T 



T=T C 



d 2 s Q 
d?T 



(36) 



T=T C 



Hence the energy density is a continuous function of T, but its first derivate has a discontinuity, for illustration see 
Fig. El right. 

What is the order of the PT in the model? To answer this question Eq. H19|) should be rewritten as: 



sh = fn + u / dvv~ a exp{-v(s H - s Q )} , 
Jv n 



(37) 



where a = —(7 + S + 2) > 1. Differentiating both sides with respect to T (the prime denotes the d/dT derivative) 
one gets: 

P 00 poo 
s 'h = f'n + u ' / dvv~ a exp [-v (s H - sq)\ +u / dvv~ a+1 (s'q-s' h ) exp [-v (s H - sq)] , (38) 



from which follows: 



G + F ■ s' Q 
1 + F ' 



(39) 
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where 

/>oo />oo 

G = f' H +u dvv- a exp[-v(s H - s Q )} , F = u dv V - a+1 exp [-v (s H - s Q )] . (40) 
Jv„ ' Jv 

It is easy to see that the transition is of the 1 st order, i.e. s'q(Tc) > s' H (Tc), provided a > 2. The 2 nd or higher 
order phase transition takes place provided s'q(Tc) = s' H (Tc) at T — Tc- This condition is satisfied when F(T) 
diverges to infinity at T — ► (Tc — 0), i.e. for T approaching Tc from below. 

One notes that the exponential factor in p(m,v) (|21J) generates the singularity sq(T) Ij28(l of the Z(s,T) function 
dl7fl . Due to this there is a possibility of a phase transition in the model. Whether it does exist and what the 
order is depends on the values of the parameters 7 and 5 in the pre-exponential power-like factor, v 1 (m — Bv) d , of 
the mass- volume bag spectrum p(m,v) 1)21(1. This resembles the case of the statistical bootstrap model discussed 
in Sec. II. The limiting temperature Tjj = 1/6 appears because of the exponentially increasing factor exp(fem), in 
the mass spectrum (J2J , but the thermodynamical behavior ((811011 at T = Tr depends crucially on the value of the 
parameter a in the pre-exponential factor, mT a (JJJ. 

What is the physical difference between a > 2 and 1 < a < 2 in the model? From Eq. I|37[l it follows that the 
volume distribution function of quark-gluon bags at T < Tc has the form 

W(v) cx v~ a exp [— v (s H — sq)] . (41) 

Consequently the average bag volume: 

poo 

dv v W(v) , (42) 



Vo 



at T — > (Tq — 0) approaches: 



v = const , for a > 2 ; (43) 

v — * 00 , for l<a<2. (44) 



Thus in the vicinity of the 1 st order PT the finite volume bags ("hadrons") dominate at T < Tc- There is a single 
infinite volume bag (QGP) at T > Tc and a mixed phase at T — Tc (see also Ref. L 2lj). For the 2 nd order and/or 
"higher order transitions" the dominant bag configurations include the large volume bags already in a hadron phase 
T <Tc, and the average bag volume increases to infinity at T — > (Tc — 0). 

The condition for the 2 nd order PT can be derived as following. The integral for the function F in Eq. I)40|l reads 

/>oo 

/ dvv- a+1 exp[~v( SH -s Q )} = (s H ~s Q y 2+a T[2-a,(s H -s Q )V ] , (45) 

where T(k, x) is the incomplete Gamma-function. Thus using Eq. <|39|l one finds at T — > (Tc — 0): 

S Q - s 'h « (s H - s Q f~ a , for a < 2 , (46) 

s'q - s' B oc - hr 1 (s H - s Q ) , for a = 2 , (47) 

and consequently: 

s' H - s' Q cx (s h - sq) 3 " 2 " for a < 2 ; (48) 

" " m ~ 3 (sg - S Q) f o , Aa s 

s H — s cx — for a = 2 . (49) 

v Sff - SQ 

Therefore for 3/2 < a < 2 the 2 nd order PT with s'^(T c ) = 00 takes place whereas for a = 3/2 the 2™ d order PT 
with the finite value of s' H (Tc) is observed. The dependence of the specific heat C = de/dT on T is shown in Fig.0J 
left. 

An infinite value of a specific heat C at T — * (Tc — 0) obtained within the model for 3/2 < a < 2 is typical for the 
2 nd order PTs jH. 

From Eq. gHl it follows that s'^(T c ) = s' Q (T c ) for a < 3/2. Using Eqs. ^1BIB5|| one finds at T -> (T c - 0): 

s'h- S '" cx - SQ ) 4 ~ 3Q . (50) 
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T c T T c T 

FIG. 4: The dependence of the specific heat on temperature for the 2 nd order PT (left) and for the 3 rd order PT (right). 
Large values of C in the hadron phase reflect the fact that large derivative den/dT at T close to Tc is needed for the 2 nd and 
3 rd order PT to reach the value of eq(Tc) without energy density discontinuity. 

Thus, for 4/3 < a < 3/2 there is a 3 rd order transition: 

sh(T c ) = SQ (T C ) , s' H (T c ) = s' Q (T c ) , s^{T c ) = s' Q (T c ) , S £(T C ) + s%{T c ) , (51) 

with s'h(Tc) — oo for 4/3 < a < 3/2 and with a finite value of s'^(Tc) for a = 4/3. The dependence of the specific 
heat C = de/dT on temperature for the 3 rd order transition is shown in Fig. 21 right. 

By calculating higher order derivatives of sh and sq with respect to T it can be shown that for (n + l)/n < a < 
n/ (n — 1) (n = 4, 5, ...) there is a n th order transition 

sh(T c ) = s Q (T c ) , s' H (T c ) = s' Q (T c ) , ... 4^ 1} (T C ) = ^^(Tc) , s%\t c ) £ 8 $>(T C ) , (52) 

with s^\Tc) — oo for (n + l)/n < a < 71/(71 — 1) and with a finite value of s^\Tc) for a = (n + l)/ri. 

The 3 rd and higher order PTs correspond to a continuous specific heat function C with its maximum at T near 
Tc (see Fig. 0J right). This maximum appears due to the fact that large values of derivative den/dT at T close to 
Tc are needed to reach the value of eq(Tc) without discontinuities of energy density and specific heat. The so called 
crossover point is usually defined as a position of this maximum. Note that in the present model the maximum of a 
specific heat is always inside the hadron phase. 



V. NON-ZERO BARYONIC NUMBER 



At a non-zero baryonic density the grand canonical partition function for the system of quark-gluon bags can be 
presented in the form 

Z(V,T^ B ) = = CX P Z ( V > T > B )> ( 53 ) 

b=~oo ^ ' 

where /is is the baryonic chemical potential (for simplicity strangeness is neglected in the initial discussion). The 
Laplace transform of Z(V,T,/ib) lEBt reads [21| 

Z{s,T,n B ) = f dVeM- sV)Z(V,T,n B ) = J— (54) 

Jo s - J{s,T,hb) 
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where 



with 



f(T,fj, B ,s) = f H (T,fj,B,s) 

/ff(T,/l B ,s) : 

p(m,v;fi B /T) \ 



dv 



v„ 



dm p(m,v; fiB /T) exp(—sv)(f)(T,m) , 



M„ + Bv 



g i exp 

OO 

S exp 



T 



bps 
T 



exp(-Vjs) 4>(T,mj) 
p(m,v;b) . 



(55) 

(56) 
(57) 



Similar to the case of ps = discussed in the previous sections one finds that the pressure is defined by the farthest- 
right singularity, s*: 



lim — \nZ(V,T,p B ) 

V— *oc V 



P{T, p B ) 

and it can be either given by the pole singularity, sh- 

s H (T,p B ) = f{s H ,T,p B ) 
or the singularity sq of the function /(s,T, pb) l|55|l itself: 



Ts*(T,p B ) = T-max{s H , s Q } , 



sq(T, p B ) 



90' 



"95 


10 




Y 


- — 


m 



5 (PbY 

9?r 4 It; 



5 

T 



(58) 



(59) 



(60) 



Note that for p B = Eq. H60|l is transformed back to sq(T) (|28|l . as (Jq(pb = 0) = ctq. The energy density and 
baryonic number density are equal to 



£(T,/i B ) = T 



2 ds*{T,p B ) 



ds*(T, p B ) 



n B (T,p B ) = T 



The second term in Eq. (|55|l . the function /q(T, ps,s), can be approximated as 



ds*(T,p B ) 
dp B 



Jq{T,s,pb) ^ u{T,p B /T) I dvv a exp[-w(s -sq(T,/x b ))] 



(61) 



(62) 



i.e., it has the same form as /q(T, s) l|29|) . At pb/T = const and T — > oo one finds sq(T, ps) oc T 3 and u(T, ps/T) oc 
yio+45 xhis is the same behavior as in the case of pb = 0, J2HJ- At a small T and ps one finds > sq, so that 
the farthest-right singularity s* equals to sjj. This pole-like singularity, sh(T, pb), of the function Z(s,T, ps) l|54l) 
should be compared with the singularity sq(T,pb) of the function Jq(T,s,Pb) (|55|l itself. The dependence of sq 
on the variables T and ps is known in an explicit form ijfJUJl. If conditions (|32[1 are satisfied it can be shown that 
with an increasing T along the lines of Pb/T — const one reaches the point Tq{pb) at which the two singularities 
collide, sq = sh, and sq(T,pb) becomes the farthest-right singularity of the function Z(s,T, ps) l|54|) at T — > oo 
and fixed pb/T. Therefore, the line of phase transitions sh{T,pb) — sq{T,pb) appears in the T — ps plane. Below 
the phase transition line one observes Sh{T,pb) > sq(T,pb) and the system is in a hadron phase. Above this line 
the singularity sq(T,pb) becomes the farthest-right singularity of the function Z(s,T, ps) (1541) and the system is in 
a QGP phase. The corresponding thermodynamical functions are given by: 



p(T,Pb) 
e(T,y B ) = T 

n B (T,p B ) = T 



5 AT " 



1627T 



4 

2 Mb 



D 



2 ds Q (T, p B ) 

dT 
ds Q (T, p B ) 
dps 



Tp B — = 3p{T,pb) + TB 



dp B 



PB T 



817T 



2 /4 



(63) 
(64) 

(65) 



The analysis similar to that in the previous section leads to the conclusion that one has the I s * order PT for a > 2 
in Eq. (JHSJ , for 3/2 < a < 2 there is the 2 nd order PT, and, in general, for (n + l)/n < a < n/(n - 1) (n = 3,4, ...) 
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there is a n th order transition. Note that sjj(T, \Ib) found from by Eq. (|59|l is only weakly dependent on a. This 
means that for a > 1 the hadron gas pressure pn = Ts^j and thus the position of the phase transition line, 

sh{T,[j, b ) = 8 Q (T,fi B ) , (66) 

in the T — fj,s plane is not affected by the contribution from the large volume bags. The main contribution to Sh 
(|59|l comes from small mass (volume) bags, i.e. from known hadrons included in fjj (|56|l . This is valid for all a > 1, 
so that the line (|66|l calculated within the model is similar for transitions of different orders. On the other hand, the 
behavior of the derivatives of Sh 1(59(1 with respect to T and/or /is near the critical line 1)66(1 . and thus the order of 
the phase transition, may crucially depend on the contributions from the quark- gluon bags with v — *■ oo. For a > 2 
one observes the 1 st order PT, but for 1 < a < 2 the 2 nd and higher order PTs are found. In this latter case the 
energy density, baryonic number density and entropy density have significant contribution from the large volume bags 
in the hadron phase near the PT line 1(66(1 . 

The actual structure of the " critical" line on the T — [Ib plane is defined by a dependence of the parameter a on 
the /is/T ratio. This dependence can not be reliably evaluated within the model and thus an external information 
is needed in order to locate the predicted "critical" line in the phase diagram. The lattice QCD calculations indicate 
that at zero hb there is rapid but smooth cross-over. Thus this suggests a choice 1 < a < 3/2 at \ib = 0, i.e. the 
transition is of the 3 rd or a higher order. Numerous models predict the strong I s * order PT at a high /is/T ratio 6], 
thus a > 2 should be selected in this domain. As a simple example in which the above conditions are satisfied one 
may consider a linear dependence, a = ao + ai/is/T, where ao = 1 + e (0 < e << 1) and a\ ss 0.5. Then the line of 
the I s * order PT at a high hb/T ends at the point hb/T rj 2, where the line of the 2 nd order PT starts. Further on 
at hb/T « 1 the lines of the 3 rd and higher order transitions follow on the "critical" line. This hypothetical "critical" 
line of the deconfinement phase transition in the T — fis plane is shown in Fig. [S] 



QGP 



1 order PT 





FIG. 5: The hypothetical phase diagram of strongly interacting matter in the T — \jlb plane within the quark-gluon bag model. 
The influence of strangeness is neglected. The line of the I s * order phase transition at a high /is/T ratio is followed by the line 
of the 2 nd order PT at an intermediate (ib/T values and by the lines of higher order PTs at a low (ab /T. 

In the case of the non-zero strange chemical potential fxs the pole singularity, sh, and the singularity sq become 
dependent on ^5. The system created in nucleus- nucleus collisions has zero net strangeness and consequently, 

ns(T,HB,^s) =T = 0. (67) 

At a small T and /is, when sh > sq, Eq. ((67(1 with s* — sh defines the strange chemical potential [is — Us (X) Ms) 
which guarantees a zero value of the net strangeness density in a hadron phase. When the singularity sq becomes 
the farthest-right singularity the requirement of zero net strangeness ((67|l with s* = sq leads to //g = fjS = hb/3- 
The functions fig(T,iJ,g) and ^(T,/j,b) = are different. Consequently, the line ((66(1 of the I s * order PT in the 
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T — fiB plane is transformed into a " strip" 2 . The phase diagram in which influence of the strangeness is taken into 
account is shown schematically in Fig. |SJ 



1 order PT 




B 



FIG. 6: The hypothetical phase diagram of strongly interacting matter in the T — [ib plane with strangeness taken into 
account. The line of the 1 st order phase transition is transformed into a "strip". For further details see caption of Fig. |S] 



The lower and upper limits of the "strip" are defined, respectively, by the conditions: 



(68) 



The T and hb values inside the "strip" correspond to the mixed HG-QGP phase for which the following conditions 
have to be satisfied: 



6- nf + (1 - £)■ nf 



mix 

n s 



5- T 



(1 - 8) ■ T 



ds 



H 



. 



(69) 
(70) 



where S is a fraction of the mixed phase volume occupied by the QGP. Eq. I|69|) reflect the Gibbs equilibrium between 
the two phases: the thermal equilibrium (equal temperatures), the mechanical equilibrium (equal pressures), the 
chemical equilibrium (equal chemical potentials). The net strangeness does not vanish in each phase separately, but 
the total net strangeness of the mixed phase (|70(l is equal to zero. This means the strangeness-anti-strangeness 
separation inside the mixed phase |2^|. The line l|66[l for the 2 nd and higher order PTs (i.e. for 1 < a < 2) remains 
unchanged as ^ becomes equal to /ig along this line. 



VI. SUMMARY AND OUTLOOK 



The phase diagram of strongly interacting matter was studied within the exactly solvable statistical model of the 
quark-gluon bags. The model predicts two phases of matter: the hadron gas at a low temperature T and baryonic 
chemical potential fiB, and the quark-gluon gas at a high T and/or fiB- The order of the phase transition is expected 
to alter with a change of the hb/T ratio. The line of the 1 st order transition at a high hb/T is followed by the line 
of the 2 nd order phase transition at an intermediate hb/T, and then by the lines of the "higher order transitions" at 
a low Hb/T. The condition of the strangeness conservation transforms the 1 st order transition line into the "strip" 
in which the strange chemical potential varies between the QGP and HG values. 



2 Within a phenomenological construction of the l 3 ' order PT this was discussed in Ref. |24| . 
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In the high and low temperature domains the approach presented here reduces to two well known and successful 
models: the hadron gas model and the bag model of QGP. Thus, one may hope the obtained results concerning 
properties of the phase transition region may reflect the basic features of nature. 

Clearly, a further development of the model is possible and required. It is necessary to perform quantitative 
calculations using known hadron states. These calculations should allow to establish a relation between the model 
parameters and the T and \ib values and hence locate the critical line in the T — fin plane. Further on one should 
investigate various possibilities to study experimentally a rich structure of the transition line predicted by the model. 
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